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Collisional dynamics of solitary matter waves of superfluid hard core bosons, consisting of dark and bright 
non-linear waves as well as supersonic periodic trains, reveals remarkable richness and coherence, with the 
phase of the condensate playing a key role. Depending upon the condensate density and their relative velocity, 
two distinctive collisional types emerge: the intuitively expected repulsive collisions due to the hard core boson 
constraint; and, also collisions in which they "pass through" one another. In addition to confirming the soliton 
status of both bright and dark solitary waves, our studies reveal a variety of multi- solitons including multiple 
families of breathers, that can be produced and precisely controlled via quantum phase engineering. 



Solitary waves and solitons are encountered in systems as 
diverse as classical water waves (il, magnetic materials 121. 
fiber-optic communication as well as Bose-Einstein con- 
densates (BEC)|4|. Rooted in nonlinearity which balances 
dispersive effects, solitons are fascinating non-linear waves 
that encode coherence underlying the system. Solitons are 
distinguished from solitary waves by their behavior in colli- 
sions. Solitons remain intact after collision, which is usually 
associated with the integrability of their non-linear dynamics. 
Therefore, discovery of solitons in a complex physical system 
suggests the existence of a hidden simplicity of the underlying 
nonlinear equations of motion. 

Intrinsically nonlinear in nature due to inter-particle inter- 
actions, the BEC systems are natural fertile ground for ex- 
ploring solitons. In dilute atomic gaseous BECs which are 
simply described in terms of the properties of the non-linear 
Schrodinger equation (NLSE), or Gross-Pitaevski equation 
(GPE), bright and dark solitons are characterized not only by 
persistent density anomalies, but also by characteristic phase 
modulations across their profiles. The quantum phase playing 
quite different roles in the bright (attractive condensates)! 5 1 
and dark (repulsive condensates), both of which have been 
experimentally realized by 161, iTTl . 

Here we report the existence of multiple species of soli- 
tons in strongly repulsive BECs as might arise in considera- 
tion of a hard core boson (HCB) system on a ID lattice |8|. 
The presence of both dark (density notch) and bright (waves 
of elevation) solitons in the same physical system is a nov- 
elty in ultra-cold atomic systems as the well known solitons in 
weakly interacting gaseous BECs are dark or bright depend- 
ing upon the repulsive or the attractive nature of the interparti- 
cle interaction.|4| Considering collisions of like pairs (bright- 
bright or dark-dark) if their relative propagation speed ex- 
ceeds a (density dependent) critical value, two colliding soli- 
ton adapt to the hard core constraint by broadening, and seem- 
ingly to pass through each other. Below the critical speed, the 
solitons bounce of each other, and their density peaks (bright- 
bright) or minima (dark-dark) repel visibly, mimicking the re- 
flection of a wave from a wall, as the hard core density con- 
straint stops them from passing through one another. In both 



types of collisions, the dynamics is controlled by the quan- 
tum phase of the composite non-linear wave and the density 
is found to remain stationary during such a pair collision. 

When the soliton speed exceeds the speed of sound, solitons 
occur as trains. When these supersonic trains collide, they 
always pass through each other. 

In contrast to the collisions involving members of the same 
soliton species ( dark-dark or bright-bright collisions), the col- 
liding solitons of different species ( dark-bright) always pass 
through each other. The special case near half-filling pro- 
vides a collision scenario resembling the soliton-antisoliton 
breather pair of the sine-Gordon system.| 1 1 In this case, the 
two species are mirror images of each other, with respect to 
the ambient surface of constant density po with respect to 
which they move, and they disappear as perturbations in the 
density at the mid-point of the collision and re-emerge un- 
scathed. Their exceptional stability is rooted in the conserva- 
tion of quantum phase jump through the collision. As we dis- 
cuss below, such inter species collisions are characterized by 
a stationary phase profile, whereas the intra species collisions 
result in stationary density profile. As in the sine-Gordon case, 
the HCB system supports breathers: oscillatory bound states 
of dark and bright solitons. In contrast to Sine-Gordon on 
the other hand, these breathers can be dissociated into dark 
and bright soliton pairs by tuning a system parameter, or as 
in the sine-Gordon case by small alterations in initial density 
or phase conditions. We show that with appropriate quantum 
phase imprinting it is possible to create a families of breathers 
where the number, and oscillation frequency, of the breathers 
can be tuned by changing the phase profile. 

The system under consideration here is the limiting case of 
the extended Bose Hubbard model in D dimensions, 
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Here, 6^ and bj are the creation and annihilation operators for 
a boson at the lattice site j, rij is the number operator, a labels 
nearest-neighbor (nn) sites , t is the nn hopping parameter. 
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U is the on-site interaction strength, and /i is the chemical 
potential. An attractive nearest-neighbor interaction V < 0, 
is introduced to soften the effect of a strong repulsive onsite 
interaction, U >> and may also describe systems with 
long range interactions such as dipolar gases.|9 | In HCB limit 
(U oo), no more than one boson can occupy a given site. 
In this case, we can represent the system by a lattice of spin- 
1/2 particles where two spin states correspond to two allowed 
boson number state |0 > and |1 >. Thus, the HCB system is 
explicitly described by. 
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where g = t — V and the spin flip operators = Sx ^ iSy 
are the annihilation and the creation operators of the corre- 
sponding bosonic Hamiltonian ,bj Sj'. Thus the order 
parameter that describes EEC wave function is i/jj = {S^'). 
In this mean-field description, the evolution equation for the 
order parameter is obtained by taking the spin-coherent state 
average of the the Heisenberg equation of motion. The spin 
coherent state \rj > at each site j can be parametrized as: 

t> +e*^ sin ^1 ;>]. With this 
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choice, the HCB system is mapped to a system of classical 
spins where the particle density pj, and the condensate 
density pj satisfy pj = PjP^, where p^ = 1 — pj is the hole 
density. In this representation, ipj = ^/pje*"^. We cast the 
equations of motion in terms of the canonical variable (j) and 
cos(6>) = (1 — 2p), where for simplicity we denote cos(^) 
equal to S are: 
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on the background particle density po implies the existence 
of a counterpart /{x^p^), superposed upon a corresponding 
hole density p^. In fact it is easy to see that /^(z,po) = 
±/^(z,Pq). For po < 1/2, the ± corresponds to bright 
and dark solitons respectively. The bright solitons have the 
unusual property of persisting at speeds upto the speed of 
sound, sharp contrast to the dark species that resemble the 
dark soliton of GPE whose amplitude goes to zero at sound 
velocity. Furthermore, unlike dark solitons, bright waves , in 
their density profile always include p = 1/2, the branch point 
of the equation ps = p{l — p). It should be noted that for 
Po > 1/2, the dark and bright solitons switch their roles. In 
view of the particle-hole duality, we will present our results 
for Po < 1/2 where bright solitons have the persistent char- 
acter noted above. Note that for v > Cg, F becomes purely 
imaginary and equation (6) then describes an infinite periodic 
train of soliton. 
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In the continuum approximation, the equations have been 
shown to support solitary waves| 8 1 riding upon a background 
density po: p{z) = po -\- f{z), with z = x — vt. 



±\/(Po -Po)2+472popScosh f - (pg - Po) ' 
Here 7 = Vl — v'^ with v the speed of the solitary wave in 



units of the speed of sound velocity c. 
is the width of the soliton, F~^ = 7 
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characteristic phase jump associated with the solitary waves 
is, A0± = v/n:2^[±sin-i +^]. 

The presence of two species of solitary waves is a di- 
rect consequence of particle-hole symmetry underlying the 
equations of motion. The existence of f{z^po) superposed 



FIG. 1 : (color online) Soliton 3D density profile (left) and its phase 
(right) projected into 2D illustrating collisions where solitons pass 
through ( top, V = .85) or repel (bottom, v = .5) each other. Here 
V = .9 and po = .45. As seen from the figure the exact soliton 
collision time can be precisely determined by monitoring the phase: 
In T-type, it becomes uniform throughout the lattice at the instant of 
maximal overlap while in R-type, the collision is signaled by a phase 
jump of TT, a phase wall, see also Figure 3B for another view of this 
phase wall. 



We now report results of a detailed study of collisions of 
identical solitary waves, and mirror image bright-dark pairs, 
moving with opposite velocities on a lattice, based on integra- 
tion of equations ^ in time, with equation ^ used to con- 
struct initial conditions. We first discuss intra species colli- 
sions of two solitary waves. Figures ([T]) present examples of 
the two distinctive types of such collisions, referred to here 
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FIG. 2: (color online) Soliton-soliton collision phase diagrams, top 
and bottom left hand figures. The dotted regions indicate parameter 
regimes where solitons engage in T class collisions (see Fig. 1), and 
appear to pass through each other; the white regions below these are 
those of the R class collisions, as in Fig. 2. The parabolic black line 
shows the sound velocity Cs as a function of p. The blue dots show 
the supersonic region, applicable only for bright soliton trains. Top 
and the bottom figures on the right respectively correspond to typical 
bright and the dark individual solitons both shown for p — A5 and 
V = .5, respectively. 



as the T, complete 'Transmission,' and R, partial 'Reflection,' 
classes respectively. These collisions are both characterized 
by stationary density , when at some time t, ^ = for all 
j. In the T-class, where solitary waves appear to pass through 
each other, the phase ^ becomes uniform across the lattice, ( 
from eq. ([5]). In the R-class, where solitary waves appear to re- 
pel each other, is uniform across the lattice except at two 
singular lattice sites, where there is an anti-node in the den- 
sity and a corresponding phase jump of tt , as encountered in 
wave reflection phenomena. The distance between the nodes 
can be tuned by changing the interaction V. Both classes dis- 
play a stationary density profile at the collision center, and 
for the R-class the particle or hole densities always attain the 
maximum density (the anti-node) permitted by the hard core 
constraint. In all of our investigations, the solitary waves were 
found to be solitons: they emerge intact after collisions. The 
numerical results presented fully validate the continuum limit 



solutions equation ([6]). The equivalent collisions for dark-dark 
pairs are not shown, but mirror those of Figs.Q , as density 
notch collisions where the anti-node maxima become actual 
nodes (zeros) in the superfluid density, and have correspond- 
ing phase jumps of tt, as would be the case for an ordinary 
node in a Schrodinger wave-function. 

We summarize the characteristics of bright and dark pair 
soliton collisions in a collision phase diagram as shown in 
Figure Q. The bright solitons repel each other ( R class) pro- 
vided their their velocity is below a threshold value. This is 
rooted in the fact that at low superfluid density, where either 
the particle or the hole population is large, the hard core con- 
straint suppresses T class collision. For dark solitons, on the 
other hand, there exists a critical velocity for T class colli- 
sions for all values of the background density. In other words, 
both bright and dark solitons encode the hard core constraint 
of the underlying Hamiltonian. The supersonic solitons are 
always in the T-class, illustrating a remarkable adaptability of 
the hard core constraint. 

In the case of half-filling, bright and dark solitons are mir- 
ror images. Then, interspecies collisions provide a conjugate 
scenario where the collision center is characterized by a sta- 
tionary phase, ^ = for all j ,rather than the stationary 
density of intraspecies collisions. The colliding solitons tem- 
porarily exactly annihilate each other in terms of density pro- 
file, but with a strongly variable phase profile, which allows 
them to re-emerge as counter propagating density defects of 
opposite sign. Reminiscent of the soliton-antisoliton pair of 
the sine-Gordon system|l |, the HCB soliton carries a topo- 
logical charge that is encoded in the the phase jump. 

As Fig. ^ shows, during a bright-dark soliton collision of 
mirror image pairs there is a smooth phase jump of approxi- 
mately 27r across the lattice even at the moment of annihila- 
tion. Thus phase imprinting of a two-7r phase jump of uniform 
gas should produce a soliton-antisoliton pair. 

The uniform density and smooth phase profile scenario as 
illustrated in ^ opens a possibility to create breathers\lO\\ 
spatially localized modes that oscillate in time and which may 
be thought of a bound state of soliton and antisoliton. We 
have found a specific phase imprinting condition that cre- 
ates such breathers. An initial condition of uniform density 
at half filling of the lattice and phase profiles of the form 
n7rtanh((x — Xi)/r) , where the integer n/2 determines the 
number of breathers in the multi-soliton profile near Xi . See 
figure (|4]). These breathers, can be dissociated into bright-dark 
soliton pairs by tuning the strength of the attractive interac- 
tion V in the Hamiltonian of Eqn 1 . Small perturbations on 
the phase profile as shown in figure ^ lead to a spectrum of 
stable oscillatory modes of sharply differing frequencies, each 
resembling a bound state of bright-dark soliton bound pair. 

The HCB constraint prohibits particles moving in one- 
dimension from passing through each other. Oddly, though, 
in the HCB systems treated here, we find solitons and soli- 
ton trains that easily pass through each other in low and high 
(supersonic) speed collisions. At low lattice fillings, the sys- 
tem also exhibits solitons that encode the fermionization char- 
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FIG. 3: (color online) Soliton density profiles (left) and phases (right) 
for T-type (a, A) and R-type (b,B) bright-bright soliton collisions, and 
a bright-dark (c,C) interspecies collision. Shown are denisties and 
phases before (red), at (black), and after (blue) the collision time, the 
blue and red curves being at times ±4 time units with respect to the 
exact collision time, as defined by the stationarity of the phase den- 
sity or phase respectively, for the intra and inter species collisions. 



acteristics of HCB, and that prohibit transmissive collision. 
Unusually, unlike the GPE solitons of gaseous BECs, both 
bright and dark solitons can be found in the same system, for a 
fixed set of parameters in the lattice Hamiltonian. Non-GPE- 
type solitary waves have been the subject of various recent 
studies. |[m[T2ll . The present studies add a new class of soli- 
tons and soliton trains relevant to strongly interacting BEG. 
The existence of solitons in mean-field equations for HGB, 
perhaps has its roots in the integrability of the underlying 
XXZ spin Hamiltonian in a magnetic field.| 13] Our numer- 
ical studies suggest the existence of a variety of novel multi- 
solitons that include double-node (or double antinode) soli- 
tons that repel, and of a broad class of breathers whose fre- 
quencies and integrations may be controlled in the laboratory 
by quantum phase engineering. 
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